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THE COMBINED EFFECT OF GRAVITY GRADIENT AND AERODYNAMIC 
TORQUES ON THE ATTITUDE CONTROL OF THE APOLLO/S- IVB 
IN A CIRCULAR ORBIT ABOUT THE. EARTH 

SUMMARY 


An analysis of the combined effect of aerodynamic and gravity gra- 
dient torques on the attitude control of the Apollo/S-IVB configuration 
during the earth-parking-orbit phase of a lunar landing mission is de- 
scribed. The parking orbit is assumed to be circular above the surface 
of a spherical earth of homogeneous mass distribution. Aerodynamic 
torques are determined by means of free-molecule impact theory, commonly 
referred to as specular reflection. 

The equations utilized in the analysis were programed on the IBM JO^k 
computer using an assumed orbit altitude of 100 nautical miles. The ef- 
fects of the disturbance torques on the limit-cycle characteristics of the 
vehicle were determined using the S-IVB attitude control system to main- 
tain the vehicle at various pitch, yaw, and roll attitudes relative to the 
local horizontal within specified deadbands. Pure pitch attitudes (pitch 
axis normal to orbit plane, zero roll and yaw) of 0, h5, and 90 degrees, 
relative to the local, horizontal were simulated using specified deadbands 
varying from ±0.5 degrees to ±5 degrees. At all commanded attitudes, the 
vehicle damped to a steady-state limit cycle on one side of the pitch 
deadband instead of traversing the entire deadband as desired. The fre- 
quencies of these limit cycles were considerably higher than those which 
would exist in the absence of disturbing torques, resulting in increased 
control system fuel consumption. 

A combined maneuver in which the vehicle was commanded to maintain 
simultaneously a pitch and yaw attitude relative to the local horizontal 
of 50 degrees each within a ±5 degree specified deadband about each axis 
was also simulated. The limit-cycles about both the pitch and yaw axes 
exhibited a similar behavior to the pure pitch command described above, 
although the amplitudes and frequencies were different in pitch from 
those in yaw. 


INTRODUCTION 


A satellite in earth orbit is acted upon by numerous environmental 
disturbance torques which affect its angular motion. These include 
torques resulting from gravity gradient, solar pressure, electric and 
magnetic fields, -meteorites, and aerodynamics. Rough estimates given by 
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Roberson in reference 1 indicate that the most significant of these 
torques for a low-altitude orbit result from gravity gradient and aero- 
dynamics . The effect of gravity gradient on the angular motions of an 
earth satellite has been analyzed by various authors such as Dool'in 
(ref. 2) . The effect of aerodynamics on the angular motions of certain 
satellite configurations at low altitudes has been studied by DeBra 
(ref. 3) . 


The present study investigates the combined effect of gravity gra- 
dient and aerodynamics on the attitude control of the Apollo/S-IVB in a 
circular orbit about the .earth. According to present plans, the Apollo/ 
S-IVB configuration will be maintained at some specified attitude rela- 
tive to the local vertical within predetermined deadbands in pitch, yaw, 
and roll by means of reaction jets during the earth orbital phase of the 
lunar landing mission. This requirement arises from the need to orient 
communications equipment on board the vehicle to earth-based flight moni- 
tors, as well as to orient an onboard sextant for star and landmark 
sightings to establish vehicle position. The S-IVB control system will 
be the primary means of attitude control and the Apollo service module 
control system may be used as a backup. 

The mathematical model employed to analyze the combined effect of 
aerodynamic and gravity gradient . torques on the attitude control of the 
vehicle, as well as results obtained using the S-IVB attitude control 
system to maintain various attitudes relative to the local vertical, is 
herein presented. 

The author wishes to acknowledge the assistance of Mr. James L. Long 
of the Computation and Data Reduction Division, NASA Manned Spacecraft 
Center, who programed and checked the digital computer mechanization of 
the equations used in this study. 


SYMBOLS 


A Length of cone removed to form cone frustum (fig. A-k ) , 

ft 

a Distance from vehicle center of mass to elemental mass 

in vehicle, ft 

C Aerodynamic force component along x-axis, lb 

C^ Normal force coefficient, dimensionless 
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Constants used to determine control torque direction 
about the x-, y-, and z-axes, respectively, (equations 60 
to 65 ) , dimensionless 

Diameter of cylinder used to represent Saturn S-IVB 
booster, ft 

Diameter of cylinder used to represent service module 
plus LEM shroud 

Differential operator, dimensionless 

Cone center of pressure distance from apex along 
x-axis, ft 

Command module center of pressure distance from nose 
along x-axis, ft 

Combined center of pressure distance from combined 
center of mass along x-s^xis, ft 

Cylinder center of pressure distance from front along 
x-axis, ft 

Cone frustum center of pressure distance from small end 
along x-axis, ft 

Large cone center of pressure location used to determine 
cone frustum center of pressure location (appendix A), ft 

Small cone center of pressure used to determine cone 
frustum center of pressure (appendix A), ft 

Service module plus lunar excursion module shroud center 
of pressure distance from front of service module along 
x-axis, ft 

Saturn S-IVB center of pressure distance from front 
along x-axis, ft 


M 


Resultant torque acting on vehicle, ft- lb 


M , M , M Control torque about x-, y-, and z-axes, respectively, 
X J y j Z J ft-lb 


Mj_, M 2 , M 3 


Magnitude of control torque about x-, y-, and z-axes, 
respectively, ft-lb 



< Ph 
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m Mass of vehicle, slugs 

N Aerodynamic normal force acting on vehicle, It 


1 

R 


r 


r cyl 


r 

m 


S 

T 


1 


t 

V 


x, y, z 


X 




z 

Q 


Unit vector perpendicular to an elemental area on surface 
of vehi cle } dimensl onle s s 

2 

Dynamic pressure, lt/ft 

Vector distance from earth center of mass to vehicle 
center of mass, ft 

Radius of cone intermediate cross-section (fig. A-3), ft 
Radius of cylinder cross-section (fig. A-5) 

Vector distance from earth center of mass to elemental 
mass in vehicle, ft 

Cross-sectional area of S-IVB, itD /k, ft 

Orthogonal matrix transformation from the x q , y Q , and. 
z^-axis system to the x, y, and z-axis system, dimen- 
sionless 

time, sec 

Velocity of vehicle relative to the atmosphere (equal to 
inertial velocity in present analysis), ft/sec 

Vehicle body axes having origin at combined center of 
mass with x-axis positive forward through command module 
nose, y- and z-axes completing right-handed system 

Cone body axes having origin at apex with x -axis 

c 

positive rearward, y - and z -axes completing right-handed 

c c 

system (fig. A- 2) 


x* , y' c > z' c Axes parallel to x c ~ and z c ~axes, respectively, with 
origin at center of elemental area on surface of cone 
(fig. A-2) 



Distance along x-axis from nose of command module to 
vehicle combined center of mass, ft 

Distance along x-axis from nose of command module to 
vehicle cc ibined center of pressure, ft 


Z cyl ( ^yli n<ier "body axes having origin at end with 


x ,-axis 
cyl 


positive rearward, y^^- and z ^.-axes completing 
right-handed system (fig. A-5) 


Inertial cartesian axis system having origin at center 
of spherical earth with y^-axis normal to orbit plane, 

x q — and z^-axes completing right-handed system (fig. 2) 


Cartesian coordinate system having origin at center of 
spherical earth with x^- axis parallel to local hori- 
zontal (positive in direction of motion) y^-axis 
normal to orbit plane and z^-axis along R (positive 
outward from center of earth) 

Combined angle of attack, or angle between x-axis and 
V, degrees 

Angle in y c z Q -plane between y c ~axis and line from 
x c ~axis to projection on y c z c ~plane of point on sur- 
face of cone (fig. A-2), degrees 

Values of p^ corresponding to points on line along 

cone surface separating portion of surface exposed to 
air flow from portion shielded from air flow, degrees 

Angle related to cone frustum corresponding to p 
for cone, degrees 

Angle related to cone frustum corresponding to p 
for cone, degrees u 

Angle related to cylinder corresponding to p for 
cone, degrees u 

• 

T t. Angle in orbit plane between z -axis and z -axis, 
O o r ’ 

degrees 



Constant orbital angular rate, or angular rate of 

x y and z -axes relative to x -, y -, and z -axes, 
r ' " r ' r o' o' o ' 

rad/sec 

Semi- apex angle of cone used to represent command 
module, degrees 

Semi-apex angle of cone frustum used to represent 
adapter: between LEM shroud and S-IVB, degrees 

Angle between V and line of intersection between cone 
surface and plane containing V and x -axis (fig. A-l), 
degrees 

Angle between V and line normal to cone surface in 

plane containing V and x -axis (fig. A-l), degrees 

c 

Euler angles defining orientation of xyz-triad relative 

to x y z -triad, degrees 
o o o 7 

Euler angles defining orientation of xyz-triad relative 
to x^y^z^-triad, degrees 

Commanded values of 0 , Y , and 0 • , respectively, 
degrees r r 

Atmospheric density, slugs/ft J 

Width (duration) of attitude control jet pulses about 
x— , y-, and z-axes, respectively, seconds 

Time between attitude control jet pulses about x-, y-, 
and z-axes, respectively, seconds 

Aerodynamic roll-angle, or angle between M and y-axis 
in yz-plane, degrees 

Resultant angular rate of vehicle in inertial space, 
rad/sec 


Subscripts 

Adapter between service module plus LEM shroud and 
S-IVB 



Aerodynamic 


Cone 

Error functions used to command control torques about 
the x-, y-, and z-axes, respectively, dimensionless 

Resultant aerodynamic force acting on vehicle, lb 

Inverse- square-law gravitational force acting on 
vehicle due to earth, lb 

Lower (negative) deadbands about the x-, y-, and z-axes 
respectively, radians 

Upper (positive) deadbands about the x-, y-, and 
z-axes, respectively, radians 

Moments of inertia about the x-, y-, and z-axes, 
respectively, slug- ft ^ 

Unit vectors along the x-, y-, and z-axes, respectively, 
dimensionless 

Unit vectors along the x c -, y Q -> and z-axes, 
respectively 

Geocentric gravitational constant of the earth, 
l.J+07656 X 10 l6 ft 5 /sec 2 

Rate gains in attitude control system, rad/rad/sec 

Attitude gains in attitude control system, rad/rad 

Length of cylinder used to represent S-IVB booster, ft 
Length of cone used to represent command module, ft 

Length of cylinder used to represent service module 
plus LEM shroud, ft 

Length of cone frustum used to represent adapter be- 
tween the service module plus LEM shroud and the 
S-IVB booster, ft 
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CM 

cyl 

f 

G 

L 

S 

sm/lem 

S-IVB 

T 

x, y, z 
x , y , 


Adapter center of pressure distance from front along 
x-axis, ft 

Command module 

Cylinder 

Cone frustum 

Gravity gradient 

Large cone 

Small cone 

Service module plus lunar excursion module shroud 
Saturn booster used for translunar orbit injection 
Total or resultant 

Components along x-, y-, and z-axes, respectively 
Components along x q -, y Q -, and z^-axes, respectively 


In addition, a ( * ) over a symbol represents a derivative with respect to 
time and a (~) over a symbol represents a vector. 


ANALYSIS 


Equations of Motion 

The mathematical model employed for this investigation was a rigid 
vehicle moving in a circular orbit about a spherical, homogeneous, non- 
rotating earth. A block diagram of the mathematical model is shown in 
figure 1. The coordinate systems used are shown in figure 2 and the 
geometry of the vehicle is shown in figure 3- The body axes of the ve- 
hicle are assumed to be principal axes and the moments of inertia of the 
vehicle are assumed to be constant. Under these conditions, the rota- 
tional equations of motion are, from reference k, 

I 1 * 


8 x * 


M + M 


G 


x. 


— Q* Q, 

y z 


- 


V 


( 1 ) 
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- ft ft (I 

X z 


X 



I 


y 


( 2 ) 


ft = Tm +M +M - ft ft (i - I )1 /l (3) 

z [_z a z G z J X y y x J / z 

The effect of attitude control jet and aerodynamic forces on the 
orbit parameters are assumed to be negligible for the length of time 
during which the vehicle remains in earth orbit, that is, the velocity 
and altitude of the vehicle are constrained to remain constant, the 
plane of the orbit being fixed in inertial space. Consequently, the 
translational equations of motion are not included in this study. 

Inertial Attitude 

The attitude of the vehicle relative to inertial space can be de- 
termined by use of Euler angles . From figure 4, the components of angular 
velocity about the body axes in terms of Euler angles and Euler angular 
rates are seen to be 

ft = 0+9 sin ijr (4) 

x 


ft 

y 


= t sin 0 + 9 cos i|r cos 0 


(5) 


ft = f cos 0-9 cos i|r sin 0 (6) 

z 

Solving equations (4), (5) , and (6) simultaneously to obtain the 
Euler angular rates. 


0 = ft - 9 sin ijf 

X 


(T) 

ilr = ft sin 0 + ft cos 0 
y r z r 


(8) 

9 = jfty cos 0 - ft z sin 0^ 

j cos \|r 

(9) 


Equations ( 7 ), (8), and (9) can then be integrated numerically with 
respect to time to obtain the attitude of the vehicle relative to inertial 
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space_, that is, the orientation of the xyz-triad relative to the 

x y z -triad, 
o^o o 


Attitude Relative to Local Vertical 

The attitude of the vehicle relative to the local vertical can he 
determined hy use of Euler angles hy superimposing a fourth angle to 
represent the rotation of the local vertical relative to inertial space, 
as shown in figure 5. As can he seen from figure 5, 


ft = (r + 0 ) sin ijr + 0 
x ' r r r r 


( 10 ) 


• • • 

ft = (r + 9 ) cos \|r cos 0 + i|f sin 0 

y r r r r r r 


ft = if r cos 0 - (r + 9 ) cos t sin 0 
z r r x ' r r r v 


(ID 

( 12 ) 


Solving equations (IQ), (ll) , and (12) simultaneously to obtain the 
relative Euler angular rates, 


K = U x~ ^ + 9 r^ Sln ^r 


(13) 


^ = ft sin 0 + ft cos 0 

r y r z r r 


( 1*0 


e r = jft y cos 0 r - ft z sin 0 r - r cos t r J j cos ( 15 ) 

Equations (l3)r (1*0* and (15) can then he integrated numerically 

with respect to time to obtain the attitude of the vehicle relative to 

the local vertical, that is, the orientation of the xyz-triad relative 

to the x y z -triad, 
r r r 


Coordinate Transformations 


The components of vehicle position and velocity along the x -, y -, 
and z Q -axes are, from figure 2, 0.0 


R = R sin r 
x 

o 

R =0 
y o 


( 16 ) 

( 17 ) 
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and. 


R 

7 

= R cos r 

( 18 ) 

O 



V 

x o 

- v cos r 

(19) 

V 

y o 

= 0 

(20) 

V 

*7 

= -V sin T 

(21) 


o 

The components of orbital angular rate along the x q -, y an< 3- 
z^-axes are, from figure 


• 

r 

- 0 

(22) 

X 

0 




= f 

(23) 

y o 

K 

= 0 

(24) 


o 


The matrix transformation from the x y z -triad, to the xyz-triad is, 
I o o o 

from figure 4, 



(cos f cos 9) 

(sin 0 sin 0 
- sin f cos 9 cps 0) 
(sin f cos 9 sin 0 
+ sin 9 cos 0) 


(sin f) 

(cos 0 cos \Jr) 

(-sin 0 cos \|r) 


(-sin 9 cos \|0 
(sin 9 sin \|r cos 0 
+ sin 0 cos 9 ) 

(cos 0 cos 9 
-sin 0 sin 9 sin i|r) 


(25) 


The components of vehicle position and velocity along the x-, y- , 
and z-axes are then 
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1 

** 

1 



1 

O 

* 

1* 

R 

y 

- 

7^ 

1 1 


R 

y o 

1 

N 

PS 

1 



1 

tsi 

0 

1 


and 


— 


— - 

V 

X 


V 

X 

0 

V 

y 

-IM 

V 

y o 

1 

1 


■ 1 

O 


( 26 ) 


( 27 ) 


The components of orbital angular rate along the xyz-axes are 



Aerodynamic Torque 


The configuration shown in figure 3 is assumed to represent the 
Apollo / S- IVB for the purpose of aerodynamic calculations . 

As shown, the command module is represented by a right circular 
cone, the service module/LEM shroud and Saturn S-IVjB by right circular 
cylinders and the S-IVB adapter by a right circular cone frustum. 


The total normal force coefficient for this configuration is 
determined by summing the normal force coefficients of the individual 
components, that is. 
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( "'E = Si + Sf / + Si + Si 

T "cm sm/lem A S-IVB 


( 29 ) 


where all coefficients are based on a common reference area and shielding 
of various components by upstream components are assumed to be negligible. 
The coefficients for the individual components are derived in Appendix A. 
The resulting expressions are listed below, each based on the cross- 

sectional area of the S-IVB, tcD /k . For the command module. 


N. 


= 2 


CM 


'D 
. D 


cos 


5^ sin 2a 


when a < 5^ 


( 30 ) 


/j) \2 

- 1 \rj 00t M I [\ + 1) sin 2a 8ln 25 i 


+ cos p. 


2 2 1 

2 cos a sin ^ s l n 2a sin 25^ sin 3^ 


2 2 2 / 2 \ 
+ sin a cos (sin 3^ + 2 

\ u J 


when a > 8, 


( 31 ) 


where, from figure 6, 


a = tan 


Tr 2 2 

•l V y + V z 
V 


( 32 ) 


and, from Appendix A, 


3 1 = cos 

u 


2 2 

^ y tan a - tan 5^ 


tan a 


( 33 ) 


For the service module/LEM shroud, 


K 


sm/lem 


32 L 2 D l .2 

T- o~ sin a 

3« d 2 


( 34 ) 



( 35 ) 


2A 


For the S-IVB adapter. 


V 2 


- (V d ) ! 


cos 62 sin 2 a when a < 8 g 


C H A - ; cot 5 2 


[1 - ( Di /d) 2 _ 


I ( B 2 + 2 ) Sin 2a Sin 2S 2 

v \ u 


+ cos P 2 


r 2 

, 2 cos a 

u*- 


2 1 

sin Sg - g sin 2 a sin 28g sin (3g 


u 


2 2 2 / 2 \* 
+ — sip a cos S 2 (sin (3g + 2 j 


when a > 8 ,. 


( 36 ) 


where, replacing (3 and 8 by (B 0 and 8 „ respectively, in equation 

(33), u 1 2 u 2 


P 2 = cos 
u 


, 2 2 
^ y tan a - tan 8 ^ 

tan a 


( 37 ) 


For the Saturn S-IVB, 


■'N, 


S-IVB 


32 

3it 


(§) sin2a 


( 38 ) 


From figure J, the resultant aerodynamic moment about the nose of 
the command module, 


M 


D 


2 _ 


a = ‘l* 31- 
NOSE 


p n 4- P 


sm/ehm ( l i + 1 sm/lem) 


+ Cjj (h^ + Lg + 1^) + Cjj ( 

A VlVB ' 


L 1 + l 2 + + 1 S _ IVB ) 

_ 

( 39 ) 



from Appendix A, 

*CM = 3 

(40) 

HlCVJ 

H 

H 


L,(D„ + 2D) 

n D JL 

'A " 3 (D x + D) 

m 

^S-IVB = 2 L 

m 


Substituting equations (4o) through (4j) into equation (39) and simpli- 
fying, 

M W 1 ^ [ (2/3> %< 

+ 1 f(l/2)C + 

{_ sm/lem 

+ ■+ 2D) p 

The resultant normal force on the vel 
\ - C N T 


+ + 

SM/LEM \ S-IVB 


+ Si 

“a “s-ivb 


+ (l/2) L C. 


% 


s-ivb 


(D 1 + D) + 0 




S-IVB 


hide 


m 


The distance of the resultant center of pressure from the nose of 
the command module is then determined by the ratio of resultant moment 
about the nose to resultant normal force, that is. 


X ='P- = XoSE / *' 


(W) 



Substituting equations (44) and (45) into equation (46) and reducing* 


X ,. P . ’ Ch 


(2/3) C + 0 + 0 + 0 

CM SM/LEM A. ^S-JVB 


+ Lr- 


+ L-, 


(l/2) C K , + °H + C N 

sm/lem a s-ivb 


+ (1/2) LC. 


X 




(D + 2D) C h (D + D) + C ]\ / 
A' 1M S-IVB [i 


% 


T 


The distance of the resultant center of pressure from the resultant 
center of mass 


1 

c 


.p. 


« X 


c.m« 


- X 


c,p. 


(48) 


The resultant aerodynamic moment about the center of mass 


M a - % * S 


c-P- 


and acts about a» axis normal to the plane containing the x-axis and the 
resultant velocity vector, as shown _n figure 6, The components about 
each of the body axes are 

M » 0 (50) 

x a 

M y = M a cos 0 b (51) 

M = M sin 0, (52) 

z a Hs ' 

a 


where, from figure 6, 
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Gravity Gradient Torque 


The torque produced by gravity gradient on a vehicle in an inverse- 
square force field is treated in Appendix B. Expressions are developed 
to determine the components of this torque about an arbitrary set of 
orthogonal body axes of an arbitrary vehicle. Since the body axes of 
the Apollo/S-IVB are assumed to be principal axes in the present analysis, 
the components of gravity gradient torque about these axes are (see 
Appendix B), 



R • R (I 
R 5 y z v z 


M 

■ y 


= 2| R R (I 
t, 5 x z v : 


G R- 


X 



f R E (I 
R 5 x y y 


V 

(5*0 


(55) 


(56) 


Attitude Control System 

The attitude control system employed in the present analysis con- 
sists of constant magnitude jets mounted on the vehicle to produce torque 
about each of the body axes . 


The jets producing control torque about each of the body axes are 
fired at constant frequency and constant width (duration) as commanded by 
the control logic for each axis. The system is assumed to be ideal, 
that is, the pulses are rectangular and no hysteresis or time delays 
exist . 


The error signals used to command control torques about the x-, y-, 
and z-axes respectively are: 

E x ■ K 1 ( fl x - E x> + K a y - K ) <57) 

c 

E y = K 3 (! V ' V + K 4 (0 r ■ s r c > < 58) 

\ " % K - q> + K 6 <*r - > 

c 


( 59 ) 
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The quantities r , T , and r in equations (57) j (58), and (59) > 
x y z 

respectively, remove the angular rate of the x^y^z^-triad relative to 
inertial space so that the error functions E , E , and E are 

x ‘y z 

analogous to a familiar inertial attitude command system rate 
feedback • 


The control torques about the x~, v-, and z-axes, respectively. 


are: 


M - C;L m, 


(6o) 


M * C 2 M 2 


(6i) 


M = CL M, 

z T 3 3 


(62) 


where 


( 1 for t seconds, then 0 for t seconds 

V ' vr 


X 


and repeat if E < H 
^ x x n 


x n 


^ -1 for seconds, then 0 for seconds, and repeat if ^(63) 


E > H 
x Xr- 


V. 

r 


0 if H < E < H 

X 1 — X - x 2 


C 2 => 


J 

'l 


1 for t seconds, then 0 for t seconds, and repeat 


yi 


if E < H 

y y x 


- 


1 for t seconds, then 0 for t seconds, and repeat if 


E > H 

y y£ 

0 if H < H j 

^ * 1 ~ *2 


.(64) 



19 


C„ - 


r ' 

1 for T seconds, then 0 for t seconds, and repeat 
z Zj. 7 

if E < H 
z z^ 

^ -1 for t z seconds, then 0 for seconds, and repeat 

if E > H 
z z„ 


> (65) 




0 if H < E < H 
Z 1 “ z " z 2 


J 


TEST PROCEDURE 


All attitude command maneuvers simulated corresponded to a Circular 
orbit at an altitude of 100 nautical miles above the surface of the 
earth. The vehicle and orbital parameters for all maneuvers were: 


A = 18.833 ft 

D = 21.66? ft 

D x = 12.833 ft 

_ , k p 

I = 8.02 X 10 slug- ft 

I = 9.75 x 10 6 slug- ft 2 

y 

I. = 9*75 X 10 6 slug- ft 2 

K = 1.4075 x 10 l6 ft 5 /sec 2 

= K,_ = 5 rad/rad/sec 

Kg = =. Kg = 1 rad/rad 

L 


= 43.625 ft 
-II.583 ft 
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Lg = 32.700 ft 

= 12.917 ft 
= 3. 489 X 10 3 ft-lb 

Mg = 3-702 X 10 3 ft-lb 

M ^ = 8. 476 X '10 3 ft-lb 

R = 2.15115 x 10 7 ft 

V = 2.5579 x 10 k ft /sec 

X = 73-10 ft 

C.Hl. 

r = 6.812468 X 10 2 deg/sec 
5^ = 33- 06 deg 

6g =18.85 deg 

p = 9.67 x 10" 13 slugs /ft 3 

T - T = t = 0.05 seconds 

X j z 

T v ' - T „ = T „ = O.95 seconds 

X I y l Z 1 

A series of pitch command maneuvers were simulated requiring the 
initial conditions listed below: 


ijf = 0 

/ = 0 



f = 0 
r 

c 
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0 =0 
r T 

C 

ft =0 

X 

ft = r = 6.812468 X lO - ^ deg/sec 

y 


r = o 


Commanded pitch attitudes (0 ) of 0, 45, and 90 degrees were 

c 

simulated. At each attitude, results were obtained for desired deadbands 
ranging from +0-5 degrees to +5 degrees. 

A combined pitch and yaw maneuver was simulated wherein the initial 
conditions listed above were applicable with the following changes: 


E 


x. 


= H = H 


-5 deg 


H x = H v “ H z = 5 de § 
x 2 y 2 z 2 


9 = 35 deg 
t ~ 35 deg 
0 =0 deg 


0 r = 35 deg 


t r = 35 deg 

0 r = 0 deg 

0 r = 30 deg 
c 

f r = 30 deg 
c 
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& = 0 
/ Vi 


ft = r = 0.03907 deg/sec 


x 


ft = r = 0 . 0^571 deg/sec 

y y 

ft z = T z = 0.03201 deg/sec 

Each of the above maneuvers was simulated for a time sufficient to allow 
the limit-cycle periods about each axis to damp to steady-state values. 


RESULTS AND DISCUSSION 


A series of maneuvers using the Saturn S-IVB attitude control system 
to maintain the Apollo/S-IVB at different pitch attitudes relative to the 
local horizontal was simulated. The vehicle was initially oriented with 
the y-axis normal to the orbit plane; consequently, the resultant aero- 
dynamic and gravity torques acted about that axis, the x- and z-axes being 
unaffected. In addition, a combined maneuver in which the vehicle was 
commanded to maintain simultaneously a pitch' and yaw attitude relative to 
the local horizontal, a situation which could be required for a navita- 
tion sighting or other operation, was also simulated. 


PITCH COMMAND 


The vehicle damped to a steady- state limit cycle on one side of the 
specified pitch deadband instead of traversing the entire deadband at 
each commanded pitch attitude. This behavior is a consequence of the uni- 
directional disturbance torques (fig. 7) acting about the y-axis tending 
to increase the pitch attitude. As an example of this behavior, consider 
the case listed in table I for a commanded pitch attitude of 0 degrees 
with a specified deadband of +0.5 degrees. As shown in figure 8, when 
the vehicle reaches the +0. 5-degree side of the deadband from within the 
deadband at a sufficiently low rate, the pitch attitude control jet will 
impart a minimum impulse to the vehicle and reverse the direction .of 
rotation. The disturbance torques are then acting in a direction oppo- 
site to the direction of rotation with sufficient magnitude to stop the 
rotation before the vehicle reaches the -0.5-degree side of the deadband 
and force it to start back toward the +0. 5-degree side of the deadband. 
This behavior continues until a steady- state condition is reached wherein 
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the vehicle rotates 0.053 degrees away from the +0. 5-degree side of the 
deadband in the direction of the -0.5-degree side, reverses direction 
and returns, the cycle repeating (and the jet firing) every 413 seconds. 
Similar results were obtained at the other commanded attitudes examined, 
as shown In table I . 

The resulting limit-cycle widths and periods were found to vary 
with specified deadband width at the commanded attitudes examined. This 
variation came about because the vehicle always remained in close prox- 
imity to one side of the specified deadband, thus subjecting the vehicle 
to varying (attitude-dependent) disturbance torques. This variation is 
seen from table I to be especially significant for a commanded pitch 
attitude of 0 degrees, the resulting limit- cycle period varying from 
413 seconds for -a +0.5 degree- specified deadband to 33 seconds for a 
+5-degree specified deadband. For this reason, a series of maneuvers 
were simulated wherein the commanded pitch attitude was 0 degrees and 
the specified deadband widths were varied from +0.5 degrees to +5 de- 
grees in increments of 0.5 degrees. Both the deadband widths and limit- 
cycle periods resulting decreased as the specified deadband widths 
increased. The results obtained are summarized in table II. 

Combined Pitch- Yaw Command 

The pitch and yaw attitudes were each commanded to 30 degrees within 
a ±5-degree deadband. The limit- cycles about both the pitch and yaw 
axes exhibited similar behavior to the pure pitch command described above. 
The pitch deadband width obtained was 0.001 degrees with a period of 
7 seconds at the +5-degree side of the deadband. The yaw deadband width 
obtained was 0.008 degrees with a period of 54 seconds at the -5-degree 
side of the deadband. The difference in limit-cycle characteristics in 
pitch and yaw was due in part to the gravity gradient and aerodynamic 
torques being in the same direction (positive) in pitch, but in opposite 
directions in yaw. Another source of .difference was the inertia cross- 
coupling between the pitch and yaw axes through the roll axis because a 
roll rate was required as a component of the orbital angular rate (rota- 
tion rate, of local horizontal in the orbit plane) . Although this cross- 
coupling was insignificant compared to the gravity gradient and aero- 
dynamic torques about the pitch axis, it dominated the motion about the 
yaw axis . The final difference arose because the control torque about 
the yaw axis was greater than the control torque about the pitch axis 
(two 150- pound-thrust jets produce torque in either direction for yaw 
control, whereas only one 150- pound- thrust jet is available for pitch 
control in each direction), while the moments of inertia in pitch and 
yaw were equal. 
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CONCLUDING- REMARKS 


Gravity gradient and aerodynamic torques acting on the Apollo/S-IVB 
in a 100- nautical- mile circular orbit about the earth are of sufficient 
magnitude to have a significant deleterious effect on the limit- cycle 
operation of the Saturn S-IVB attitude control system about both the pitch 
and yaw axes of the vehicle. At the commanded attitudes considered in the 
present analysis, the torques forced the vehicle to damp to a steady-state 
limit-cycle on one side of the pitch and yaw specified deadbands instead 
of traversing the entire deadbands', resulting in frequencies considerably 
higher than the no-disturbance-torque condition. 

At commanded attitudes near zero, the frequency of the steady-state 
limit-cycles about the pitch or yaw axis can be considerably reduced by 
decreasing the specified deadbands, resulting in a decrease in attitude 
control fuel consumption. At larger commanded attitudes, the effect of 
varying the specified deadbands is less pronounced. 
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TABLE I.- EFFECT OF VARIATION IN COMMANDED PITCH ATTITUDE ON THE 
PITCH-AXIS LIMIT-CYCLE CHARACTERISTICS OF THE APOLLO/S- IVB 
IN A 100- NAUTICAL-MILE CIRCULAR EARTH ORBIT* 


Commanded 
pitch attitude 
(degrees) 

Specified 
deadband width 
(degrees) 

Resulting 
deadband width 
(degrees) 

Resulting 
limit cycle 
period 
(sec) 

o 

±5 

.003 

33 

45 

±5 

.001 

5 

90 

±5 

.002 

14 

0 

±0.5 

• 053 

413 

45 

±0.5 

.001 

6 

90 

±0.5 

.001 

11 


*S-IVB attitude control system with minimum pulse duration of 
50ms and 150- pound- thrust jets used. 
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TABLE II.- EFFECT OF VARIATION IN SPECIFIED PITCH DEADBAND ON THE 
PITCH- AXIS LIMIT- CYCLE CHARACTERISTICS OF THE APOLLO/S- IVB 
IN A 100- NAUTICAL-MILE CIRCULAR EARTH ORBIT* 


Commanded 
pitch attitude 
(degrees) 

Specified 
deadband width 
(degrees) 

Resulting 
deadband width 
(degrees) 

Resulting 
limit cycle 
period 
(sec) 

0 

±0.5 

0.053 

413 

0 

±1.0 

0.027 

199 

0 

±1.5 

O.OI7 

131 

0 

±2.0 

0.013 

98 

0 

±2.5 

0.011 

78 

0 

±3.0 

0.008 

64 

0 

±3-5 

0.007 

55 

0 

±4.0 

0.006 

48 

0 

±4.5 

0.006 

43 

0 

±5.0 

0.003 

33 


*S-IVB attitude control system with minimum pulse duration of 
50ms and 150- pound- thrust jets used. 
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Figure 1.- Block diagram of mathematical model 


















Figure 4.- Euler angles defining attitude of 
vehicle relative to inertial space. 




Figure 5 — Euler angles defining attitude of vehicle 
relative to local vertical. 
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Figure 6 Orientation of resultant velocity and aerodynamic 
torgue vectors relative to vehicle body axes . 




Roll axis displacement from local horizontal, deg 
Figure 7.- Disturbance torques aeting on Apollo/S-IVB in 

1 m-rmn-M r»«l —mi 1 ^ nirnnlar a"hrM vf- -Kha dqt41o 
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Figure 8.- Comparison of desired and achieved pitch-axis phase plane trajectory of 
Apollo/S-IVB in a 100- nautical- mile circular orbit about the earth. 
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APPENDIX A 

THEORETICAL STATIC AERODYNAMIC FORCE COEFFICIENTS AND 
CENTERS OF PRESSURE FOR CERTAIN SYMMETRICAL BODIES 
OF REVOLUTION IN FREE-MOLECULE FLOW . 

•RIGHT CIRCULAR CONE 


Force Coefficients 

Consider a body having a perfectly smooth surface to be moving in 
a highly rarefied gas, that is, a gas in which the frequency of colli- 
sions between individual molecules is negligibly small. According to 
reference 5? if a sufficient number of molecules are present to determine 
the macroscopic properties of the gas (pressure, temperature, density, 
et cetera), the condition is described as free-molecule flow. Near the 
surface of the body the two streams of incident and emergent molecules 
experience little interaction. In this case, the air molecules are said 
to reflect specularly; that is, the molecules are assumed to be smooth, 
perfectly elastic spheres which will rebound from a smooth surface at 
the same angle at which they strike. The behavior is illustrated in 
figure A-l for a stream of air molecules impinging oh a surface element 
of a right circular cone. 


z 

c 



Figure A-l 
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In order to deflect the stream of incident molecules through the angle 
2e shown in figure A-l, the surface element ds must impart a velocity 
component 


Vp = 2 J V j cos T] = 2V cos T] 


(Al) 


to the. stream perpendicular to the surface element. If the velocity 
change is assumed to occur in unit time, the force exerted on the sur- 
face element due to the momentum change of the stream of molecules 


where 


dF = MV 
P 


(A2) 


M 


pV ds 
H P 


(A3) 


Substituting equations (Al) and (A3) into equation (A2), 


dP 


it-pV 2 


cos T]ds 




But p and V can be expressed in terms of dynamic pressure 

<1 = | pv 2 (A5) 

Substituting equation (A5) into equation (A4) , 

dF = 8q cos 2 Tlds (a6) 

For a right circular cone, the position of a surface element Is shown in 

figure A- 2. 

z 

c 



Figure A-2 
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In order to determine the aerodynamic force acting on the element, the 
angle T| can he expressed a? follows . From figure A-2, a unit vector 
perpendicular tp the element of surface area 


P = P . i +P.j + P , k 
x' c y' "c z' c 
c v c c 


(A7) 


and the velocity vector 


V = V , i +V, j + V , k 
x' c y c z' c 
c ■ C c 


(A8) 


The components of P along the body axes are 


P - P sin 8. 
x' 1 


P , = -P cos cox 


P , = -P cox 8 sin p n 


J 


and the components of V along the same axes are 

V . = V cos a 

•v- * 


V . * 0 


V z , = v Sin a 
Q 


J 


(A9) 


(A10) 


Since Tj is the angle between P and V, it can be obtained from the 
dot product 


V • P = VP cos TJ = V , P, + V , P, + V , P , 

1 x* x' y y z' z 1 
c c “c “ q c c 


(All) 


Substituting equations (A9) and (A10) into equation (All), 


VP 


cos 71 * VP cos a sin 8^ - VP cos 8^ sin p.^ sin a. 
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or 


cos T) = cos a, sin 


cox 8^ sin p^ sin a 


(A12) 


Substituting equation (A12) into equation (a6), the resultant force 
acting on the surface element is 


2 

dF = 8q (cos a sin 8^ - cos 8^ sin p^ sin a) ds (A13) 

The resultant force can be divided into a normal force dN- along the 

z -axis and an axial force dC along the x -axis, 
c c 


dN = -cLF cos 5^ sin P^ 

(Al4) 

dC = dF sin 5^ 

(A15) 


If the body possesses homogeneous mass distribution, the center of mass 
is located on the axis of symmetry, or the x^axis of figure A-l. Since 

only the moment produced by the resultant aerodynamic force is of interest 
in the present analysis and since the axial force acts through the center 
of mass, the axial force equation, equation (A15), will not be developed 
further . 


The element of area shown in figure A-l can be seen to be 

rap n dx 

dS = 

cos 8^ 

Substituting equation (Al6) into equation (Al4), 

dN = dFr sin pdp, dx fdS 
1 c 


(ai6) 


(A 17 ) 


When a < 8^, the entire surface of the cone is exposed to the flow. 

However, when a > 8^, a portion of the surface is shielded from the 

flow. According to impact theory, the pressure in the shielded region 
is zero and thus contributes nothing to the aerodynamic force . The 
boundary line separating the exposed region from the shielded region is 
the locus of points along the surface at which the velocity is tangent 
to the surface. Thus, at these points, the velocity vector V is per- 
pendicular to the area p and cos 7] =0. 



To determine P^ at points along the boundary, equation (A12) can be 
written 


cos 71 = cos a. sin 8, - sin a cos -8 n sin p, = 

■u 1 1 1 


0 


(A18) 


u 


or rearranging and solving equation (Al8) for P^ 


u 


, tan 5 

3, .= sin" - 

1 tan a 

u 


(A19) 


The resultant normal force acting on the body can be found by integrating 
equation (A17) over the exposed surface area, that is. 


x 3- 

u 1 


N = - 


0 


J dPr sin p^P-jdx^i dS 


(A20) 


Substituting equation (A13) into equation (A20), 

X c P 1 

p u p u 2 

N = -8q / / (cos a sin 5^ - sin a cos 5^ sin- p^) sin p^d0.jdx 

0 "I 

2 (A21) 

In order to integrate equation (A2l), it is convenient to define a 
normal force coefficient per unit length of the body. 


C* 





where 



JL = _J_ o 

qS qJtr^ 


or, differentiating equation (A23), 



dN 

2 

q^ 


(A22) 


(A23) 


(A24) 



Substituting equation (A2^) into equation (A22) , 


4l 


l 


dN 


N„ _ 2 d(x /D, ) 

C qotr, ' c' 1 


(A25) 


Differentiating equation (A2l) with respect to x, 

K 


dW 

dx 


n u 2 

- 8qr j (cos a sin - sin a cos 8^ sin 3-^) sin P^P-P^ 


-at 

2 


(A26) 


since r is a function of x c only. 
But 


dJ[ _ 1_ dl 
dx " d 1 


(A27) 


Thus, the right side of equation (A26) can be multiplied by D^ and the 
resulting quantity substituted into equation (A25), yielding 


N, 


P 1 

-8qrD 1 n u 

/ (cos a sin 5 - sin a cos 5 sin p,) sin 3..dp, (A28) 

aatr. '-i-at 1 . ± 


G qatr 1 u -rt 


Substituting Dg = 2r^ and expanding equation (A28), 


C« 


= z l6.r_ r 
H C ' * r l i. 


U (cos^a sin^ 5^-2 sin a cos a sin 6^ cos 5^ sin 3^ 


2 2 2 

+ sin a cos sin 3^_) sin P^^^^ 


(A29) 


Since the cone angle 5^ is constant, straightforward integration of 
equation (A29) holding a constant during the integration yields 



k2 


c 'k„ ■ I f, {I (s u + sin 2a sln 28 i + oos p i 

C 1 I u| 


2 2 
2 cos a sin &.. 


1 2 2 2 2 
2 sin 2a sin 26^ sin 0^ + r? sin a cos 6^ (sin 0^ + 2) 

u ^ u 


(A30) 


For the case in which a < 5^, the entire surface of the cone is exposed 
to the flow. In this case the upper limit of integration is 0 1 

and equation (A30) reduces to 


1 2 
u 


k r 

C' = — sin 2a sin 28., 
N c r 1 1 


(A31) 


When a > 8^, equation (A19) must he utilized to evaluate equation (A30) 
It can he seen from equation (AI 9 ) that 

tan 2 a 


tan 2 6, 


cos 0 


1 tan a 

u 


(A32) 


Substituting equation (A32) into equation (A30) and rearranging, 


n c = ! ^ cos2 8 i sin 2a ^ P i u + tan 6 i + 5 


cos 0 


u 


cot a tan 5, 


+ 2 tan a 


(A33) 


It should he remembered that equations (A3l) and (A33) were derived on 
the assumption that centrifugal forces resulting from flow over a curved 
surface are negligible and that expansion flow regions contribute nothing 
to the aerodynamic force. 

The total normal force coefficient can be obtained from equation (A22), 



k-3 



(A3*0 


When a < 5^_, substitution of equation (A3l) into equation ( A3^) yields 


Cjj =4 sin 2a sin 25 1 



(A35) 


But, as shown in figure A- 3, 


r 



x 

c 



(A36) 


z 

c 



Figure A- 3 


Multiplying and dividing equation (A36) by 



(A37) 



Also from figure A- 3, 


Thus 



2 tan 8^ 


(A38) 


f- = 2 # tan 8 (A3 9) 

1 1 


Substituting equation (A39) into equation (A35).> replacing sin 28.. by 
2 sin 8^ cos 8^, 



1 6 


sin 


& 1 sin 2a 




(A4o) 


■ Integrating equation ( A4o) and evaluating at the upper and lower limits. 


=8 

C 


2 

sin 8^ sin 2a 


(Akl) 


But, from equation (A 38 ), 

L 1 1 cos 8 1 
2 sin 8.^ 

Substituting equation (A42) into equation ( Akl ) , 
2 

CLg. ' = 2 cos 8, sin 2a when a < 8 n 


(A42) 


(a43) 


When a > 8^, substitution of equation (A39) into equation (A30) followed 
by substitution of the resulting expression into equation (A34) yields, 
after integration along x q , 
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% • ; cot 6 


i ^l( p i u + f) sin 2a sln 2S i + cos \ 


2 2 
2 cos a sin 6.. 


1 2 2 2 / 2 \ 
2 sin 2a sin 25^ sin p^ + j sin a cos 6^ I sin p^ + 2j 

^ u D \ u / 


when a > 5., 


(Akk) 


It should be remembered that C as evaluated from equation (a 43) or 

C 2/ 

(Akk) is based on the maximum cross-sectional area of the cone, 


Center of Pressure 

In order to determine the aerodynamic moment on a right circular cone 
from the normal force coefficient, it is necessary to know the location 
of the center of pressure. This can be determined from the expression 


1 


C 


moment M 

normal force “ N 


(All-5) 


The normal force on a unit length of the cone is, from equation (A25) 


dW 



( a!+6) 


and the resulting moment about the cone vertex 

/x_ 


dM = x C 1 d _ 
c N c \D. 


•Y 


if " 1 


Integrating equation (AVf) over the length of the cone, 



(A47) 


(a48) 


0 
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Substituting equations (A3l) and (A3 9 ) into equation ( A48) , 

h 

M = 16 D 1 qitr ^ 2 sin 2<x sin^ (^j ' d^ 


(A4 9 ) 


Integrating equation (A4 9 ), 


„ 16 2 , . 2 R / n 

M = — D 1 qmr.^ sin 2 a sm 5.^ — 


(A50) 


The resultant normal force on the cone 


n - c Hc q jtr 1 


(A51) 


Substituting equations (A4l), (A50 ) , and (A5l) into equation (a 45) and 
simplifying. 


"c = 3 L l 


(A52) 


where 1 ^ is measured from the cone vertex. 


RIGHT CIRCULAR CORE FRUSTUM 


Normal Force Coefficient 

The normal force coefficient for a right circular cone frustum can 
be obtained from the expressions for the normal force coefficients of a 
right circular cone developed above. This is accomplished by subtracting 
the coefficient of the small cone shown in figure A-4 from the coeffi- 
cient of the large cone . When a < 82 , the coefficient of the small 
cone is seen from equation (A43) and figure A-4 to be 

= 2 cos 2 sin 2a (A53) 

s 

where is based on the area jtD^^/4. 



kl 




Figure A* 4 


Similarly, the coefficient fpr the large cone 


Cjj = 2 cos 5 2 sin 2a 
L 


(A5^) 


where is based on the eirea 


Before subtracting C w from C, T to determine the coefficient of the 

% A L 2/, 

frustum, both mu?t be based on the same area. If the area rtD/4 is 

chosen, must be multiplied by the ratio (jtD^A-)/ (nlF/k) = D^/d 2 , 


\\ 2 

C N S " 2 \T/ COS 5 2 sln 2a 


(A55) 


Then the frustum coefficient 




(A56) 
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Substituting equations (A 54 ) and (A 55 ) into equation (A56 ) , 

2 \ 


c w =2 

f 


1- V 

\ D 1 


cos , 5 ^ sin 2a a < 


(A57) 


When ux > 8^? the above procedure can be applied to equation (A 44 )., 

replacing 8^ . and P^ by S 2 and P 2 respectively in equation (A 44 ) 
to obtain. 


u 


u 


\ = * COt 5 2 


D. 


1 - 


_1 

D , 


\ ($2 + § \ sin 2<x sin 28 2 
' V u / 


+ COS p r 


2 2 1 

2 cos a sin ^2 ~ 2 s ^ n 2a S ^ n 28 2 S ^ n ^2 


u L_ 


2 2 2 2 
+ *7 sin a cos S 0 (sin 0 ( 

5 2 


t +2 | a 


u 


> B„ 


(A 58 ) 


•where P 2 is obtained from equation (A 32 ) by replacing 8^ and p^ 
u u 


by 8 p and P p respectively. It should be remembered that C as 
2 u 2 W f 

determined from equation (A 57 ) or (A58) is based on the area jtD / 4 . 


Center of Pressure 

The center of pressure of the frustum can be determined using 
equation (A52) with the appropriate subscripts for the small and large 
cones of figure 4 . The center of pressure of the large cone 


1 L 5 (A + L ^) 

and the center of pressure of the small cone 


(A59^ 


1 S = 3 A 


(A60) 



ij-9 


The center of pressure distance from the cone vertex is determined 
hy the ratio of large-cone moment minus small-cone moment to large- cone 
normal force minus small- cone normal force. Symbolically, 


A + 1 


f 


n l- n s 


(a6i) 


But 


and 




T\T n JtD 

W L = C N L q T" 
rtD 2 

W S ^ 

S J 




JtD 




jtD n 


M s N s 1 s V q 4 1 s 

o 


j 


(A62) 


(A63) 


Substituting equations (A62) and (A 63 ) into equation (A6l) and simpli- 
fying. 


A + 1 


f 



(a64) 


■ 65 ) 


Substituting equations (A59) and (A6o) into equation (A64) and simpli- 
fying, 
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From equation (A4-3), replacing 5. ; , "by 5g, 


C„ = 2 cos - 5 p sin 2a 

L ■ * ; . 


C w = 2 cos 8 p sin 2a 
S j 


(A66) 


Substituting equation (A66) into equation (A 65 ) and simplifying 


2 V 2 (A + Lj) - D x 2 A 
f ^ (D 2 - D 1 2 ) 


- A 


(A 67 ) 


From figure A- 4, using a property of similar triangles, 


A A + L 5 


D, 


D 


(A68) 


or, solving for A, 


A - ^ 

D-D, 


(A 69 ) 


Substituting equation (A 69 ) into equation (A 67 ) and rearranging, 


L 3 (D 1 + 2D) 
X f " 3 (Dl[ + ,D) 


(A70) 


RIGHT CIRCULAR CYLINDER 


Normal Force Coefficient 

The normal force coefficient per unit length of a right circular 
cylinder can be obtained from equations (A19) and (A 30 ) developed above 
for a right circular cone. The procedure is simply to equate the cone 
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angle 8^ in equations (A19) and (A30) to zero and determine the re- 
sulting expressions. From equation (AI9), with 8^ = 0, 

* <, 

0 3 = 0 (A71) 


Substituting equation (A7l) into equation (A30) , with 8^ = 0 
r = r^ = r ^ l° r the cylinder. 


and 


C* 



32 . 

sin 


a 


(A72) 


Using the nomenclature shown in figure A-3, 




y cyl 


the normal force coefficient for the complete cylinder 


N 


cyl 


*/ 


J cyl 

W P , 


X 


0 


cyl 

N cyl V D cyl y 


(A73) 


Substituting equation (A72) into equation (A73) and integrating, 

cyl cyl 


(A74) 
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where. 

cyl 


is based pn the area 



Center of Pressure 


The elemental torque due to the normal force acting on a unit 
length of the cylinder 


dM = X . C* d 
cyl cyl . EL 



qrrtr 


cyl 


(A75) 


and the resultant torque 


M = D . 
cyl cyl 



X £Ji] d | fc£l 


N C \ D cyl 


D I qitr cyl 
cyl/ J 


(A 76 ) 


Substituting equation (A 72 ) into equation (A76) and performing the 
indicated integration. 


w 16 cyl 2 .2 

M cyl = > 1’ a ' oy l slna 


(A7T) 


The resultant normal force acting on the cylinder 

„ 2 
ET , = C. T qitr n 

cyl N C yl cyl 

Substituting equation (A 7 ^) into equation (A78) 


(A 78 ) 


N 


32 cyl 2 . 2 

= •#— qirr sxn a 


cyl 3 « D 


cyl 


cyl 


(A79) 


The center of pressure distance from the front of the cylinder 

M . 

_ cyl 


cyl Ef 


cyl 


(A 80 ) 


Substituting equations (A 77 ) and (A 79 ) into equation (A80) and simpli- 
fying, 

hyl * 2 L cyl < A81) 
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APPENDIX B 


GRAVITY GRADIENT TORQJE ACTING ON A VEHICLE 
IN AN INVERSE SQUARE FORCE FIELD 


The force due to gravity acting on an element of mass dm in the 
vehicle shown in figure B-l can he written vectorially as: 


Kdmr 

af = - -y4 (Bl) 

r 



But, hy vector addition (see fig. B-l), 

r m =R+a (B2) 

where 

R = R x i + R y j + R z k (B3) 
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a = xi + yj + zk 


(b4) 


Substituting equations (B3) and. (b 4) into equation (B2) , 

AAA 

r = (R + x)i + (R + y) j + (R + z)k 
m v x v y . ' z 

The incremental torque on the vehicle is, from reference 6, 

dM = a xd F 

tx Gr 


(B5) 


(b6) 


Substituting equations (Bl), (B4) , and (B5) into equation (b 6) and using 
matrix notation 


,rr -Kdm 

* 3 

r 
m 


A 

(i) 

(R + x) 


J 

(y) 

(R y + y) 


A 

k 

(z) 

(R + z) 
v z 


(B7) 


Expanding equation (B7) and simplifying. 




-Kdm 


m 


1 (yE z - sB y ) + j (*®. - + £ K - y\) 


(B8) 


In order to integrate equation (B8) over the mass of the vehicle, 
it is necessary to express the quantity l/r3 as a function of R and a. 

An approximate expression of this form can be obtained in the following 


manner . 


The dot product 


r .r = r = (R + a) . (R + a) 
m m m 


= R.R+a.a+2R.a (B9) 

Substituting equations (B3) and (b 4) into equation (B9) and expanding, 

2 


m 


= R 2 1 + (a/R) 2 + (2/R 2 ) (xR + yR + zR ) (BIO) 

L 1 x y z 
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Assume now that (a/R) « 1 and can be ignored in equation (BIO) . 

Thus, 


R ) (xR + yR + zR ) 
' v x y z' 


r m 2 ~ « 2 [l + ( 2 / 

Raising equation (Bll) to the -3/2 power, 

-3 ~ -^[l + (2 /r 2 ) (xR x + yR. 
r m 


) 


(Bll) 


+ zR ) 

y z 


,-3 

2 


(B12) 


Expanding equation (B12) by the binomial theorem, retaining only the 
first term of the expansion. 


■ 1 ? ~ -4 1 - \ (xB + JE + ZE ) 

r 5 - E 3 L R 2 X y 

m 


(B 13 ) 


Combining equations (B 8 ) and (B13), the components of dM along the 
x, y, and z-axes, respectively, are: 


dM 


-Kdm 


X G R 5 


1 - - 7 : (xR + yR + zR ) 
.2 ' x J y z' 


R 


yR - zR 
z y 


(Bl4) 


dM = 

y G . 


-Kdm 


1-4- (xR„ + yR + zR ) 


R 


zR - xR 
x z 


(B15) 


dM 


-Kdm 
S G R 5 


1 - -4r (xR + yR + zR ) 
.2 ' x ' tr ^ 


R 


y 


xR - yR 

y x 


Expanding equation (Bl4), 


dM 




-Kdm 

E 3 


tv 


zR - (xyR R^ - xzR R 
R 


y t,2 x z x y 


(Bl6) 


+ y^R R - yzR ^ + yzR ^ - z^R R ) 
yz^y J z y z'' 


(BIT) 
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The total torque about the x-axis due to gravity gradient for any 
orientation of the vehicle relative to the attracting body can be found 
by integrating equation (EL7) over the mass of the vehicle with R , 

and R held constant during the integration. Integrating equation (B17) 
z 

term by term. 


M 

X G 


-K 


R 


, < R / ydm - R 
3 \ Z J y 


-R R / xzdm + R R 
X y J y Z 




+ R / yzdm - R R 
z ,/ y z 



(B18) 


Since the origin of the xyz-axis system is located at the center of mass 
of the vehicle (fig- B-l), it must necessarily follow that: 





(B19) 


to be consistent with the definition of center of mass. Substituting 
equation (B19) into equation (Bl8), then adding and subtracting 

,5, „ „ r 2. 


(30n k/v. 


x dm to the resulting expression yields: 



R R 
x y 


/ 


xzdm 


+ R R 

y z 


dm - R R 

y z 


fo 


z 2 ) 


dm 


But 



I 

I 

f 


xydm = 
xzdm = 
yzdm = 





T y z J 


(B20) 


(B21) 
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J (y 2 + z 2 ) dm s 

^ (x 2 + z 2 ) dm = X^ 
J"' (x 2 + y 2 ) dm = I z 




Substituting equations (B2i) and (B22) into equation (B20), 


m =2| rri -rri +(r 2 -r 2 )i +RR(I - I 
x G R 5 I X z xy X y xz z y yz y z' z y 


3 


(B22) 


(B25) 


For the speciai case in which x, y, and z are principal axes of the 
vehicle, equation (B23) reduces to: 


M =^RR (I - I) 
X G R 5 y z Vz y' 


(B24) 


Applying the above procedure to the y-axis, it can be shown that: 


3K 


[° 


M =^|(R -R )l +RR (I -I) 
y G r 5 I x z xz x z v x z' 


+RRI -RRI 
x y yz y z xz 


■] 


(B25) 


or, if the x, y,- z-axes are principal axes, equation (B25) reduces to: 


M «~R8 (I - I) 
y G r 5 x z x z' 


(B26) 


Similarly, it can be shown that 


M .*[■(■■ 

z e s 5 L y 


R ) I +RR (I - I ) + R R I -RRI 
x xy x y y x yzxz xzyz 


(B27) 



58 


or, if the x, y, z-axes are principal axes, equation (B27) reduces to: 



3K 

R 5 


R R 
x y 


(I 


y 




(B28) 


It should be remembered that the results, obtained above are based 
on the assumption that the gravitational force generated by the attract- 
ing body on an element of mass in the vehicle is represented by equa- 
tion (Bl) . The results are further limited by the validity of neglecting 
higher order terms in the binomial expansion used to obtain equation (B13) • 



